Traffic phenomena in biology: from molecular motors to organisms* 
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Traffic-like collective movements are observed at almost all levels of biological systems. Molecular 
motor proteins like, for example, kinesin and dynein, which are the vehicles of almost all intra- 
cellular transport in eukayotic cells, sometimes encounter traffic jam that manifests as a disease 
of the organism. Similarly, traffic jam of coUagenase MMP-1, which moves on the collagen fibrils 
of the extracellular matrix of vertebrates, has also been observed in recent experiments. Traffic- 
like movements of social insects like ants and termites on trails are, perhaps, more familiar in our 
everyday life. Experimental, theoretical and computational investigations in the last few years 
have led to a deeper understanding of the generic or common physical principles involved in these 
phenomena. In particular, some of the methods of non-equilibrium statistical mechanics, pioneered 
almost a hundred years ago by Einstein, Langevin and others, turned out to be powerful theoretical 
tools for quantitative analysis of models of these traffic-like collective phenomena as these systems 
are intrinsically far from equilibrium. In this review we critically examine the current status of our 
understanding, expose the limitations of the existing methods, mention open challenging questions 
and speculate on the possible future directions of research in this interdisciplinary area where physics 
meets not only chemistry and biology but also (nano-)technology. 

I. INTRODUCTION 

Motility is the hallmark of life. What distinguishes a traffic-like movement from all other forms of movements is 
that motile elements move on "tracks" or "trails". However, in sharp contrast to vehicular traffic, the tracks and 
trails, which are the biological analogs of roads, can have nontrivial dependence on time during the typical travel time 
of the motile elements. What makes biological traffic even more unusual is that in many cases the motile elements 
themselves not only create the tracks but also modify their lengths as well as shape and, in some extreme cases, even 
leave behind a trail of destruction by wiping out the track as they move forward. 

We are mainly interested in the general principles and common trends seen in the mathematical modeling of 
collective traffic-like movements at different levels of biological organization [l| . Although the choice of the physical 
examples and modelling strategies are biased by our own works and experiences, we put these in a broader perspective 
by relating these with works of other research groups. We begin at the lowest level, starting with intracellular 
biomolecular motor traffic on filamentary rails. Then we present brief summaries of recent works on the traffic of 
molecular motors along the collagen fibrils in the extra-cellular matrix and those on transport of micron-size cargo 
by uni-cellular micro-organisms. We end our review by discussing the collective traffic- like terrestrial movements of 
social insects, particularly, ants, on their trails. 

II. THEORETICAL APPROACHES 

In recent years many individual-based models of biological traffic have been formulated in discretized space. While 
in some models the dynamics of the system has been formulated in terms of differential equations assuming continuous 
time, in many others the dynamics has been implemented in terms of "update rules" in discrete time steps in the 
spirit of cellular automata and lattice gas models [1, |3, 0| ■ Most of these models are extensions of a class of particle- 
hopping models which were earlier successfully used in the context of vehicular traffic [s*]. The asymmetric simple 
exclusion process (ASEP) ^] is one of the simplest particle- hopping models. In the ASEP particles can hop (with some 
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probability or rate) from one lattice site to a neighbouring one, but only if the target site is not already occupied by 
another particle. "Simple Exclusion" thus refers to the absence of multiply occupied sites. Generically, it is assumed 
that the motion is "asymmetric" such that the particles have a preferred direction of motion. 

For such driven diffusive systems the boundary conditions turn out to be crucial. If periodic boundary conditions 
are imposed, i.e., the sites 1 and L are made nearest-neibours of each other, all the sites are treated on the same 
footing. If the boundaries are open, then a particle can enter from a reservoir and occupy the leftmost site (j — 1), 
with probability a, if this site is empty. In this system a particle that occupies the rightmost site (j = L) can exit 
with probability j3. The ASEP has been studied extensively in recent years and is now well understood [1, Q . Both 
parallel and random-sequential updating schemes have been studied extensively in the literature. 

The average number of motile elements that arrive at (or depart from) a fixed detector site on the track per unit 
time interval is called the flux. One of the most important transport properties is the relation between the flux and 
the density of the motile elements; a graphical representation of this relation is usually referred to as the fundamental 
diagram. If the motile elements interact mutually only via their steric repulsion their average speed v would decrease 
with increasing density because of the hindrance caused by each on the following elements. On the other hand, for a 
given density c, the flux J is given by J = cw(c), where w(c) is the corresponding average speed. At sufficiently low 
density, the motile elements are well separated from each other and, consequently, v{c) is practically independent of 
c. Therefore, J is approximately proportional to c if c is very small. However, at higher densities the increase of J 
with c becomes slower. At high densities, the sharp decrease of v with c leads to a decrease, rather than increase, 
of J with increasing c. Naturally, the fundamental diagram of such a system is expected to exhibit a maxium at an 
intermediate value of the density. 

III. INTRA-CELLULAR TRAFFIC OF CYTOSKELETAL MOLECULAR MOTORS 

Intracellular transport is carried by molecular motors which are proteins that can directly convert the chemical 
energy into mechanical energy required for their movement along filaments constituting what is known as the cy- 
toskeleton [1, Q. Three superfamilities of these motors are kinesin, dynein and myosin; majority of these motors are 
two-headed. Most of the kinesins and dyneins are like "porters" in the sense that these move over long distances 
carrying cargo along the filamentary tracks without getting completely detached; such motors are called processive. 
On the other hand, the conventional myosins and a few unconventional ones are nonprocessive; they are like "rowers" . 

These cytoskeleton-based molecular motors play crucially important biological functions, e.g., in axonal transport 
in neurons. The mechano-chemistry of single cytoskeletal motors and the mechanism of their motility have been 
investigated both experimentally and theoretically for quite some time [10. ITTI. [l2j. 

Often a single microtubule (MT) is used simultaneously by many motors and, in such circumstances, the inter- 
motor interactions cannot be ignored. In this article we shall focus mostly on the effects of mutual interactions of 
these motors on their collective spatio-temporal organisation and the biomedical implications of such organisations. 
Fundamental understanding of these collective physical phenomena may also expose the causes of motor-related 
diseases (e.g., Alzheimer's disease) [H, [13, [IE Ell thereby helping, possibly, also in their control and cure. The bio- 
molecular motors have opened up a new frontier of applied research — "bio-nanotechnology" . A clear understanding 
of the mechanisms of these natural nano-machines will give us clue as to the possible design principles that can be 
utilized to synthesize artificial nanomachines. 

Derenyi and collaborators [l^. [Tsl developed one-dimensional models of interacting Brownian motors. They mod- 
elled each motor as a rigid rod and formulated the dynamics through Langevin equations for each such rod assuming 
the validity of the overdamped limit; the mutual interactions of the rods were incorporated through the mutual 
exclusion. 

The model considered by Aghababaie et al.[l^ is not based on TASEP, but its dynamics is a combination of 
Brownian ratchet and update rules in discrete time steps. In this model the filamentary track is discretized and the 
motors are represented by field-driven particles in the spirit of the particle-hopping models. The hopping probabilities 
of the particles are obtained from the instantaneous form of the local time-dependent potential. No site can accomodate 
more than one particle at a time. Each time step consists of either an attempt of a particle to hop to a neighbouring 
site or an attempt that can result in switching of the potential from flat to sawtooth form or vice- versa. Both forward 
and backward movement of the particles are possible. However, neither attachment of new particles nor complete 
detachment of existing particles were allowed. 

The fundamental diagram of the model (Toj . computed imposing periodic boundary conditions, is very similar to 
those of TASEP. This observation indicates that further simplification of the model proposed in ref. is possible 
to develope a minimal model for interacting molecular motors. Indeed, the detailed Brownian ratchet mechanism, 
which leads to a noisy forward-directed movement of the field-driven particles in the model of Aghababaie et al. 
[Ill, is replaced in some of the more recent theoretical models [10, [Sj, lH, [H, [H [H [H, [13, [Hi by a TASEP-like 
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probabilistic forward hopping of self-driven particles In these simplied versions, none of the particles is allowed 
to hop backward and the forward hopping probability is assumed to capture most of the effects of biochemical cycle 
of the enzymatic activity of the motor. The explicit dynamics of the model is essentially an extension of that of the 
asymmetric simple exclusion processes (ASEP) (^,1131 that includes, in addition, Langmuir-like kinetics of adsorption 
and desorption of the motors. 




FIG. 1: A schematic description of the TASEP-like model introduced in ref. ^23] for molecular motor traffic. Just as in TASEP, 
the motors are allowed to hop forward, with probability q. In addition, the motors can also get "attached" to an empty lattice 
site, with probability A, and "detached" from an occupied site, with probability D from any site except the end points; the 
rate of attachment at the entry point on the left is a while that at the exit point on the right is /3. 

In the model of Parmeggiani et al. (23|, the molecular motors are represented by particles whereas the sites for the 
binding of the motors with the cytoskeletal tracks (e.g., microtubules) are represented by a one-dimensional discrete 
lattice. Just as in TASEP, the motors are allowed to hop forward, with probability q, provided the site in front is 
empty. However, unlike TASEP, the particles can also get "attached" to an empty lattice site, with probability A, and 
"detached" from an occupied site, with probability D (see figH]) from any site except the end points. The state of the 
system was updated in a random-sequential manner. Carrying out Monte-Carlo simulations of the model, applying 
open boundary conditions, Parmeggiani et al.[23| demonstrated a novel phase where low and high density regimes, 
separated from each other by domain walls, coexist [H, H^. Using a mean- field theory (MET), they interpreted this 
spatial organization as traffic jam of molecular motors. 

A cylindrical geometry of the model system was considered by Lipowsky, Klumpp and collaborators [13, [HI, [l^] to 
mimic the microtubule tracks in typical tubular neurons. The microtubule filament was assumed to form the axis of 
the cylinder whereas the free space surrounding the axis was assumed to consist of Nch channels each of which was 
discretized in the spirit of lattice gas models. They studied concentration profiles and the current of free motors as 
well as those bound to the filament by imposing a few different types of boundary conditions. This model enables one 
to incorporate the effects of exchange of populations between two groups, namely, motors bound to the axial filament 
and motors which move diffusively in the cylinder. They have also compared the results of these investigations with 
the corresponding results obtained in a different geometry where the filaments spread out radially from a central 
point. 

A novel feature of the model of Klein et al. [3l| (see Fig. [J) is that the lattice site at the tip of a filament is 
removed with a probability W per unit time provided it is occupied by a motor; the motor remains attached to the 
newly exposed tip of the filament with probability p (or remains bound with the removed site with probability 1 — p). 
Thus, p may be taken as a measure of the processivity of the motors. This model clearly demonstrated a dynamic 
accumulation of the motors at the tip of the filament arising from the processivity. 



A. Traffic of interacting single-lieaded motors KIFIA 

The models of intracellular traffic described so far are essentially extensions of the asymmetric simple exclusion 
processes (ASEP) % ^ that includes Lan gmuir-like kinetics of adsorption and desorption of the motors. In reality, 
a motor protein is an enzyme whose mechanical movement is loosely coupled with its biochemica l cy cle. In a recent 
work |32|, we have considered specifically the single-headed kinesin motor, KIFlA l33l. \3^. [35l . [37| ; the movement 
of a single KIFlA motor was modelled earlier with a Brownian ratchet mechanism [3^ |39[. In contrast to the earlier 
models [21,,^,^,^] of molecular motor traffic, which take into account only the mutual interactions of the motors, 
our model explicitly incorporates also the Brownian ratchet mechanism of individual KIFlA motors, including its 
biochemical cycle that involves adenosine triphosphate (ATP) hydrolysis. 

The ASEP-like models successfully explain the occurrence of shocks. But since most of the bio-chemistry is cap- 
tured in these models through a single effective hopping rate, it is difficult to make direct quantitative comparison 
with experimental data which depend on such chemical processes. In contrast, the model we proposed in ref. (3^ 
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FIG. 2; A schematic description of the model suggested by Klein et al. [sH for motor induced depolymerization of cytoskeletal 
filaments. 

incorporates the essential steps in the biochemical processes of KIFlA as well as their mutual interactions and involves 
parameters that have one-to-one correspondence with experimentally controllable quantities. Thus, in contrast to the 
earlier ASEP-like models, each of the self-driven particles, which represent the individual motors KIFlA, can be in 
two possible internal states labelled by the indices 1 and 2. In other words, each of the lattice sites can be in one of 
three possible allowed states (Fig.|3]): empty (denoted by 0), occupied by a kinesin in state 1, or occupied by a kinesin 
in state 2. 
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FIG. 3: A 3-state model for molecular motors moving along a MT. denotes an empty site, 1 is K or KT and 2 is KD. 
Transition from 1 to 2, which corresponds to hydrolysis, occurs within a cell whereas movement to the forward or backward 
cell occurs only when motor is in state 2. At the minus and plus ends the probabilities are different from those in the bulk. 

Good estimates for the parameters of the model could be extracted by analyzing the empirical data [s^l . Assuming 
that each time step of updating corresponds to 1 ms of real time, we performed simulations upto 1 minute. In the 
limit of low density of the motors we have computed, for example, the mean speed of the kinesins, the diffusion 
constant and mean duration of the movement of a kinesin on a microtubule from simulations of our model; these are 
in excellent quantitative agreement with the corresponding empirical data from single molecule experiments. 

Using this model we have also calculated the flux of the motors in the mean field approximation imposing periodic 
boundary conditions. Although the system with periodic boundary conditions is fictitious, the results provide good 
estimates of the density and flux in the corresponding system with open boundary conditions. 

In contrast to the phase diagrams in the a — /3-plane reported by earlier investigators [2l], [2^, [2^ , we have drawn 
the phase diagram of our model in the oja — uJh plane by carrying out extensive computer simulations for realistic 
parameter values of the model with open boundary conditions. The phase diagram shows the strong influence of 
hydrolysis on the spatial distribution of the motors along the MT. In particular, the position of the immobile shock 
depends on the concentration of the motors as well as that of ATP; the shock moves towards the minus end of the MT 
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with the increase of the concentration of kinesin or ATP or both. The formation of the shock has been estabUshed by 
our direct experimental evidence; our findings on the domain wall arc in qualitative agreement with the corresponding 
experimental observations [s^ ]. 

This work has been discussed in more detailed in our separate article (4Q] in this proceedings. 



IV. INTRA-CELLULAR TRAFFIC OF NUCLEOTIDE-BASED MOTORS 



Helicases and polymerases are the two classes of nucleotide-based motors that have been the main focus of experi- 
mental investigations. In this section, we discuss only the motion of the ribosome along the m-RNA track. Historically, 
this problem is one of the first where TASP-like model was successfully applied to a biological system. 

In a living cell ribosomes translate the genetic information 'stored' in the messenger-RNA (mRNA) into a program 
for the synthesis of a protein. mRNA is a long (linear) molecule made up of a sequence of triplets of nucleotides; 
each triplet is called a codon (see Fig. |4]). The genetic information is encoded in the sequence of codons. A ribosome, 
that first gets attached to the mRNA chain, "reads" the codons as it translocates along the mRNA chain, recruits 
the corresponding amino acids and assembles these amino acids in the sequence thereby synthesizing the protein for 
which the "construction plan" was stored in the mRNA. Once the synthesis is completed, the ribosome gets detached 
from the mRNA. Thus, the process of "translation" of genetic information stored in mRNA consists of three steps: (i) 
initiation: attachment of a ribosome at the "start" end of the mRNA, (ii) elongation: of the polypeptide (protein) as 
the ribosome moves along the mRNA, and (iii) termination: ribosome gets detached from the mRNA when it reaches 
the "stop" codon. 




FIG. 4: The process of biopolymerization: Ribosomes attach to mRNA and read the construction plan for a biopolymer which 
is stored in the genetic code formed by the sequence of codons. None of the codons can be read by more than one ribosome 
simult aneously. 

In order to model the traffic of ribosomes on a m-RNA track, let us denote each of the successive codons by the 
successive sites of a one-dimensional lattice where the first and the last sites correspond to the start and stop codons. 
The ribosomes are much bigger (20-30 times) than the codons. Therefore, neighbouring ribosomes attached to the 
same mRNA can not read the same information or overtake each other. In other words, any given site on the lattice 
may be covered by a single ribosome or none. Let us represent each ribosome by a rigid rod of length ij.. If the rod 
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FIG. 5: Schematic representation of the one-dimensional burnt bridge model of MMP-1 dynamics proposed in ref.'s^. 

representing the ribosome has its left edge attached to the i-th site of the lattice, it is allowed to move to the right by 
one lattice spacing, i.e., its left edge moves to the site i + 1 provided the site i + ir is empty. In the special case ir = 1 
this model reduced to the TASEP. Although the model was originally proposed in the late sixties [4l|, significant 
progress in its analytical treatment for the general case of arbitrary could be made only three decades later; even 
the effects of quenched disorder has also been considered in the recent literature ^H, |4^, 13, H^, . 

As mentioned above, a ribosom is much bigger than a base triplet. However, modifying the ASEP by taking into 
account particles that occupy more than one lattice site does not change the structure of phase diagram [4l| . 

V. EXTRACELLULAR TRANSPORT: COLLAGEN-BASED MOTORS 

The extracellular matrix (ECM) [47[ of vertebrates is rich in collagen. Monomers of collagen form a triple-helical 
structure which self-assemble into a tightly packed periodic organization of fibrils. Cells residing in tissues can secret 
matrix metalloproteases (MMPs), a special type of enzymes that are capable of degrading macromolecular constituents 
of the ECM. The most notable among these enzymes is MMP-1 that is known to degrade collagen. The collagen fibril 
contains cleavage sites which are spaced at regular intervals of 300 nm. The coUagenase MMP-1 cleaves all the three 
a chains of the collagen monomer at a single site. 

Breakdown of the ECM forms an essential step in several biological processes like, for example, embryonic devel- 
opment, tissue remodelling, etc. Malfunctioning of MMP-1 has been associated with wide range of diseases |48{ . 
Therefore, an understanding of the MMP-1 traffic on collagen fibrils can provide deeper insight into the mechanism 
of its operation which, in turn, may give some clue as to the strategies of control and cure of diseases caused by the 
inappropriate functions of these enzymes. 

Saffarian et al. [49| used a technique of two-photon excitation fluorescence correlation spectroscopy to measure the 
correlation function corresponding to the MMP-1 moving along the collagen fibrils. The measured correlation function 
strongly indicated that the motion of the MMP-1 was not purely diffusive, but a combination of diffusion and drift. 
In other words, the "digestion" of a collagen fibril occurs when a MMP-1 executes a biased diffusion processively 
(i.e., without detachment) along the fibril. They also demonstrated that inactivation of the enzyme eliminates the 
bias but the diffusion remains practically unaffected. They claimed that the energy required for the active motor-like 
transport of the MMP-1 comes from the proteolysis (i.e., degradation) of the collagen fibrils. 

There is a close relation between the traffic of MMP-1 on coUagens and the "burnt-bridge model" introduced by 
Mai et al. In the burnt bridge model (see Fig. [5]), a "particle" performs a random walk on a semi-infinite 

one- dimensional lattice that extends from the origin to -|-oo. Each site of the lattice is connected to the two nearest 
neighbour sites by links; a fraction c of these links are called "bridges" and these are prone to be burnt by the random 
walker. A bridge is burnt, with probability p, if the random walker either crosses it from left to right or attempts 
to cross if from right to left [13, [5l|. In either case, if the bridge is actually burnt, the walker stays on the right of 
the burnt bridge and cannot cross it any more. The hindrance against leftward motion, that is created by the burnt 
bridges, is responsible for the overall rightward drift of the random walker. Mai et al.^50] studied the dependence 
of the average drift velocity v on the parameters p and c by computer simulation. They also derived approximate 
analytical forms of these dependences in the two limits p ^ 1 and p ~ 1 using a continuum approximation. 

Saffarian et al. [49*1 also carried out computer simulations of a two-dimensional model of the MMP-1 dynamics on 
collagen fibrils which is essentially a two-dimensional generalization of the burnt-bridge model. By comparing the 
results of their simulations with their experimental observations, Saffarian et al. they concluded that the observed 
biased diffusion of the MMP-1 on collagen fibrils can be described quite well by a Brownian ratchet mechanism [3R[39|. 
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VI. CELLULAR TRAFFIC 

A Mycoplasma mobile (MB) bacterium is an uni-ccUular organism. Each of the pear-shaped cells of this bacterium 
is about 700 nm long and has a diameter of about 250 nm at the widest section. Each bacterium can move fast on glass 
or plastic surfaces using a gliding mechanism. In a recent experiment [s^ narrow linear channels were constructed 
on lithographic substrates. The channels were typically 500 nm wide and 800 nm deep. Note that each channel was 
approximately twice as wide as the width of a single MB cell. The channels were so deep that none of the individual 
MB cells was able to climb up the tall walls of the channels and continued moving along the bottom edge of the walls 
of the cannels. In the absence of direct contact interaction with other bacteria, each individual MB cell was observed 
to glide, without changing direction, at an average speed of a few microns per second. 

When two MB cells made a contact approaching each other from opposite directions within the same channel, one 
of the two cells gave way and moved to the adjacent lane. However, in a majority of the cases, two cells approaching 
each other from the opposite directions simply passed by as if nothing had happened; this is because of the fact that 
the width of the channel is roughly twice that of the individual MB cell. Moreover, when two cells moving in the 
same direction within a channel collided with each other, the faster cell moved to the adjacent lane after the collision. 

Hiratsuka et al.fl^] attached micron-sized beads on the MB cells using biochemical technique and demonstrated 
that the average speed of each MB cell remained practically unaffected by the load it was carrying. In contrast to the 
nonliving motile elements discussed in all the preceedings sections, the cells are the functional units of life. Therefore, 
the MB cells have the potential for use in applied research and technology as "micro-transporters". More recently, 
the unicellular biflagellated algae Chlamydomonas reinhardtii (CR), which are known to be phototactic swimmers, 
have been shown to be even better candidates as "micro-transporter" as these can attain average speeds that is about 
two orders of magnitude higher than what was possible with MB cells (5^. However, to our knowledge, the effects of 
mutual interactions of the CR cells on their average speed at higher densities has not been investigated. 

VII. TRAFFIC IN SOCIAL INSECT COLONIES: ANTS AND TERMITES 

From now onwards, we shall study traffic of multi-cellular organisms, particularly, ants which are social insects. 
The ability of the social insect colonies to function without a leader has attracted the attention of experts from 
various disciplines [H, [H, [H, [13, H^, [H, [13, HH ■ Insights gained from the modeling of the colonies of such insects are 
finding important applications in computer science (useful optimization and control algorithms) [6^ . communication 
engineering [6^ j_artificial "swarm intelligence" [G^I and micro-robotics [65] as well as in task partitioning, decentralized 
manufacturing [H, [g^, [H, [H, [t^, [tJ] and management [ill . the collective terrestrial movements of ants have close 
similarities with the other traffic-like phenomena considered here. When observed from a sufficiently long distance 
the movement of ants on trails resemble the vehicular traffic observed from a low flying aircraft. 

Ants communicate with each other by dropping a chemical (generically called pheromone) on the substrate as they 
move forward [zl, [tI, [T^I . Although we cannot smell it, the trail pheromone sticks to the substrate long enough for 
the other following sniffing ants to pick up its smell and follow the trail, [t^. Both the continuum model developed 
by Ranch et al. 76] and the C A model introduced by Watmough and Edelstein-Keshet [77| were intended to address 
the question of formation of the ant-trail networks by foraging ants. Couzin and Franks [7a ] developed an individual 
based model that not only addressed the question of self-organized lane formation but also elucidated the variation 
of the flux of the ants. 

In the recent years, we have developed discrete models [t^, [s^, [Ml, [12] that are not intended to address the question 
of the emergence of the ant-trail [83j] , but focus on the traffic of ants on a trail which has already been formed. We have 
developed models of both unidirectional and bidirectional ant-traffic by generalizing the totally asymmetric simple 
exclusion process (TASEP) [1, [83,[8^ with parallel dynamics by taking into account the effect of the pheromone. 

In our model of uni-directional ant-traffic the ants move according to a rule which is essentially an extension of 
the TASEP dynamics. In addition, a second field is introduced which models the presence or absence of pheromones 
(see Fig. [6]). The hopping probability of the ants is now modified by the presence of pheromones. It is larger if a 
pheromone is present at the destination site. Furthermore, the dynamics of the pheromones has to be specified. They 
are created by ants and free pheromones evaporate with probability / per unit time. Assuming periodic boundary 
conditions, the state of the system is updated at each time step in two stages. In stage I ants are allowed to move 
while in stage II the pheromones are allowed to evaporate. In each stage the stochastic dynamical rules are applied 
in parallel to all ants and pheromones, respectively. 

One interesting phenomenon observed in the simulations is coarsening. At intermediate time usually several non- 
compact clusters are formed. However, the velocity of a cluster depends on the distance to the next cluster ahead. 
Obviously, the probability that the pheromone created by the last ant of the previous cluster survives decreases with 
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FIG. 6: Schematic representation of typical configurations of the uni-directional ant-traffic model. The symbols • indicate the 
presence of pheromone. This figure also illustrates the update procedure. Top: Configuration at time t, i.e. before stage I of 
the update. The non-vanishing probabilities of forward movement of the ants are also shown explicitly. Middle: Configuration 
after one possible realisation of stage I. Two ants have moved compared to the top part of the figure. The open circle with 
dashed boundary indicates the location where pheromone will be dropped by the corresponding ant at stage II of the update 
scheme. Also indicated are the existing pheromones that may evaporate in stage II of the updating, together with the average 
rate of evaporation. Bottom: Configuration after one possible realization of stage II. Two drops of pheromones have evaporated 
and pheromones have been dropped/reinforced at the current locations of the ants. 



increasing distance. Therefore clusters with a small headway move faster than those with a large headway. This 
induces a coarsening process such that after long times only one non-compact cluster survives. 

A similar coarsening phenomenon has been observed also in the bus- route model [H, [s^]. In fact, the close relation 
between our model of uni-directional traffic on ant-trails and the bus- route model has been pointed out earlier [88| . 
In the bus route model, each bus stop can accomodate at most one bus at a time; the passengers arrive at the bus 
stops randomly at an average rate A and each bus, which normally moves from one stop to the next at an average 
rate Q, slows down to q, to pick up waiting passengers [sgI. [STj. 

In vehicular traffic, usually, the average speed of the vehicles decreases monotonically with increasing density 
because the inter-vehicle interactions tend to hinder each other's movements. In contrast, in our models of ant-traffrc 
the average speed of the ants varies non-monotonically with their density over a wide range of small values of / 
because of the coupling of their dynamics with that of the pheromone. This uncommon variation of the average speed 
gives rise to the unusual dependence of the flux on the density of the ants in our models of ant-traffic. Furthermore, 
the flux does not exhibit the particle-hole symmetry which is a characteristic of the TAEP. Details of the models and 
results on both uni-directional and bi-directional traffic of ants on trails are given in the article of John et al. in this 
proceedings [89| . The experimental data reported in the pioneering experimental work of Burd and collaborators [9^ 
were too scattered to test our theoretical predictions. However, more accurate recent data [oll, HI] establish both the 
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non-monotonic variation of the average speed with density as well as the formation of cluster by the ants. 



VIII. SUMMARY AND CONCLUSION 

In this article we have reviewed our current understanding of traffic- like collective phenomena in living systems, 
starting from the smallest level of intra-cellular bio-molecular motor transport and ending at the level of the traffic of 
social insects like, for example, ants. We have restricted our attention to those theoretical works where, in the spirit 
of particle-hopping models of vehicular traffic, the language of cellular automata or extensions of TASEP has been 
used. The success of this modelling strategy has opened up a new horizon in traffic science and, we hope, we have 
provided a glimpse of an exciting frontier of interdisciplinary research. 
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